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Abstract 

The slow-to-start models are a classical cellular automata model in simulating vehicle traffic. 
However, to our knowledge, the slow-to-start effect has not considered in modeling pedestrian 
dynamic. We verify the similar behavior between pedestrian and vehicle, and propose an new 
lattice gas (LG) model called the slow reaction (SR) model to describe the pedestrian's delayed 
reaction in single-file movement. We simulate and reproduce the Seyfried's field experiments 
at the research centre Jiilich, and use its empirical data to validate our SR model. We compare 
the SR model with the standard LG model. We test different probability of slow reaction p s in 
SR model and found the simulation data of p s = 0.3 fit the empirical data best. The RMS error 
of mean velocity of SR model is smaller than that of standard LG model. In the range of p s = 
0.1 ~ 0.3, our fundamental diagram between velocity and density by simulation coincides with 
field experiments. The distribution of individual velocity in fundamental diagram in SR model 
agrees with the empirical data better than that of standard LG model. In addition, we observe 
the stop-and-go waves and phase separation in pedestrian flow by simulation. We reproduced 
the phenomena of uneven distribution of interspaces by SR model while the standard LG model 
did not implement. The SR model can reproduce the evolution of spatio-temporal structures of 
pedestrian flow with higher fidelity to Seyfried's experiments than the standard LG model. 

Keywords: Cellular automata model, Slow reaction model, Slow-to-start model, Field 
experiment, Single-file pedestrian flow. 



1. Introduction 

In recent years, the modeling of pedestrian flows has become one of the most exciting topics. 
It has attracted considerable attention from the physical science, traffic engineering, computer 
engineering, or even social psychology science yj]. The lattice gas (LG) model, as a cellular 
automata (CA) model, is a simple but effective method to simulate the movement of a large 
number of pedestrians. With very low cost in terms of computer simulation time, it is able 
to simulate the empirical results qualitatively, even quantitatively with acceptable error. It can 
reproduce several typical self-organization phenomena observed in real pedestrian flow steadily, 



"Corresponding author. E-mail address: fangjun06@mails.tsinghua.edu.cn 
Preprint submitted to Physica A 



March 8, 2013 



such as phase transition, scaling behavior, faster-is-slower effect, stop-and-go wave and lane 
formation. 

Muramatsu introduced LG model to simulate the counter flow in the channel of subway in 
earlier times Q. Jiang extended the standard LG model through setting the maximum velocity 
more than one cell per step fl. Kuang and Fukamachi studied the binary mixture flow with 
two different velocities in their own paper Some models introduced the surrounding 

environment perception and extended pedestrians' visual field to longer distance J6(-[8|. Yu and 
Yang considered the right-moving reference of Chinese pedestrians in their own paper [0 [Toll . 
Li and Fukamachi introduced the behavior of position exchange and sidle into the LG model 
1 1 ll ll2 L Weng divided pedestrian movement into three basic behaviors, i.e., move, avoid and 
swirl lOh . The models extension of Matsui |14j] and Song II 1 311 be long to finer discretization 
or multi-grid and Song's model incorporated the social force model II 1611 into LG model. Some 
modified models dealt with evacuation dynamic in the state of emergency, such as crawl behavior 
in fire disaster 01711 and indoor evacuation without visibility 11811 . 

Calibrating and validating pedestrian models are essential to flow prediction and pedestrian 
facilities management. However, so far there have been only limited attempt about these work. 
Seyfried and his colleagues at research centre Jiilich conducted many solid field experiments 
since the mid-2000s 019W22I1 . which were performed within the Hermes projects 02311 . The tra- 
jectories of all pedestrians were generated automatically from video footage [24], which can be 
used for the calibration. 

The slow-to-start models are a classical CA model in simulating motor traffic, which were 
proposed in the 1990s and researched broadly after then. The related papers have been quoted 
extensively, such as I25U27I1 . However, to our knowledge, the slow-to-start effect has never 
applied to pedestrian dynamic modeling. Although the difference of velocity and size between 
the vehicle and pedestrian is one or two orders of magnitude, we think some similar behavior 
exists in pedestrian movement. We call this effect the slow reaction of pedestrian. From 01911 and 
videos in web site 02811 . we observe that for pedestrian flow at high density, the individual velocity 
is small and the concentration of pedestrians is sometimes reduced. In this case, pedestrian 
often make delayed reaction on the movement of the pedestrian in front and cannot follow the 
front pedestrian immediately. Similar phenomena also occur in highways traffic B25I427I1 despite 
within a smaller spatio-temporal structure. In addition, The quickness degree of reaction varies 
from pedestrian to pedestrian. Not all pedestrians within a group react to movement of front 
pedestrians simultaneously. Some individuals leave much larger than average gaps in front and 
the distribution of individual interspace is uneven. This phenomena is also reported in 02911 and 
is reproduced by social force model in 113 Oil simulating the same scenario with our work. 

Based on above reasons, we propose an new LG model called the slow reaction (SR) model to 
describe the pedestrian's delayed reaction. We try to reproduce Seyfried's experiment of circular 
passageway about single-file pedestrian flow 01911 and use its empirical data to validate our SR 
model. In addition, we compare our simulation data with standard LG model to justify our SR 
model of pedestrian flow. 



2. Models 

Our model is designed to reproduce the single-file movement in an approximate elliptical 
circle, which was the experimental set-up in Seyfried's field experiments in d (see Fig. Q}. 
The experimental area is transformed to a rectangular grid with periodic boundary as our virtual 
passageway (see Fig. The virtual passageway is composed of 43 square cells with size of 
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0.4m*0.4m, corresponding to total length l p = 17.2m in physical unit (a litter shorter than 17.3m 
in lfl9lo . The cell size of our model is the same as the standard LG model jTf] and can hold one 
virtual pedestrian at each time step. At the center of passageway, there is a measured section 
with length l m = 5cell, corresponding to the /,„ = 2m in physical unit, just like in fl9tl . All pedes- 
trians move in single-file pattern and are forbidden to overpass. Therefore, the width of virtual 
passageway is insignificant in simulation. As doing in 111 911 . we calculate the one-dimensional 
density in the longitudinal distance from the measured section. 



2.9m 




2m —31 
Measured Section 



Fig. 1. Seyfried's experimental set-up in |19| 
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Fig. 2. Virtual passageway of our simulation with periodic boundary. 



2.1. Standard LG model 

In LG model, each one moves one cell at most per time step. The positions of all pedestrians 
are updated in parallel style. The states of pedestrians in standard LG model can be divided into 
two classes according to the relative positions of successive pedestrians: the free state and the 
block state. For Free state, the gap between the target pedestrian and its front pedestrian is not 
less than one cell. We adopt the free moving velocity Vf ree = 1.24m/s in Jl9ll as our calibration 
parameter. Therefore, one simulation time step corresponds to 0.4/1. 24=0. 32s in physical unit. 
The pedestrian gets into the block state when there is no gap between it and its front one. 

2.2. SR model 

For pedestrian flow at high density, the individual velocity is small and the concentration 
of pedestrians is sometimes reduced. In this case, pedestrian often make delayed reaction on 
the movement of the pedestrian in front and cannot follow the front pedestrian immediately. In 
addition, the quickness degree of reaction varies from pedestrian to pedestrian. Not all pedes- 
trians within a group react to movement of front pedestrians simultaneously. Some individuals 
leave much larger than average gaps in front at the border of and the distribution of individual 
interspace is uneven. 
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We propose an new LG model called the slow reaction (SR) model to simulate this effect. 
Our SR model is based on the spatial distribution of adjacent pedestrians. The sphere of action is 
within one cell distance, i.e. 0.4m. Leaving just one cell from the front one, the pedestrian will 
get into slow reaction state. If the gap is more than one cell or zero, the target pedestrian will 
get into the free state or block state as the definition of 2. 1 subsection. We introduce a parameter 
called the probability of slow reaction p s to describe the quickness degree of reaction to front 
pedestrians. If someone gets into the slow reaction state, he will not follow the front moving 
pedestrian timely. The pedestrian could move forward with probability p s or remain at the same 
position with probability 1 - p s . Especially, the SR model becomes the standard LG model at 
p s - 1. The position update of SR model are calculated according to the following equations: 

( x n r d'l = 0, block state 
x" +l - I x" + p s , d" — 1, slow reaction state (1) 
\ + 1, other, free state 

where x" and v" denote the position and velocity of pedestrian ; at time n respectively, d" denotes 
the gap between two successive pedestrians (not including the length of pedestrian) and p s is the 
probability of slow reaction within the limit of < p s < 1 . 



3. Simulations 



In Seyfried's experiment [19], the testers should have been distributed uniformly before start- 
ing to move. However, due to the spatial discretization of LG model, it is difficult to distribute 
all virtual testers uniformly in the grid. It is impossible when the total number of pedestrians 
N > 22 under the length of passageway l p = 43cell. We adopt a simple method of initialization 
at the beginning of simulation. We put all testers one by one with no gap from the left boundary 
of virtual passageway. We define the circle from the header entering the upper edge of measured 
section (17cell — >18cell) to the last pedestrian going out of the lower edge (22cell — >23cell). To 
reduce the effect of nonuniform distribution, we only use the data of velocity and density from 
the 50th circle to 100th circle to obtain the fundamental diagram. 

We firstly investigate the relation between velocity and density. Through the computer simu- 
lation, we can obtain the global velocity and density in the whole passageway as well as the local 
velocity and density in the measured section. However, as no global velocity data was reported 
in 11911 . our global data cannot be compared. Therefore, we only compare the local velocity and 
density with Jl9Tl . 

Similar to the statistic method of Jl9ll . we calculate the individual velocity v, = l,„/(t° l " - 1'"), 
where & and f ut is the entrance and exit times of measured section for pedestrian i. The mean 

N N 

velocity of one cycle v = 4? 2 v* . The momentary density p{t) — 2 0, (f) jl m , where 
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The cycle density p is the mean value of the momentary density p(t) during one cycle. Be- 
cause Il9ll only reported the results with six different total number of pedestrians N = 1, 15, 
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Fig. 3. The Root-Mean-Square error e as the function of p s . 



20, 25, 30 and 34, we simulate five scenes (N — 15, 20, 25, 30 and 34) to compare with the 
empirical data in 11911 . We set p s from 0.1 to 1 with increasing step Ap s — 0.1 at each simulation 
implement. We calculate the value of density and velocity in each cycle from the 50th to 100th 
circle and get the mean value. The Root-Mean-Square (RMS) error e is calculated as follows: 



j 5 



where k represents the five scenes, v* is the mean velocity of cycles from our simulation and v™ fl " 
is the mean velocity from the manual procedure in IU9II . The RMS error e as the function of p s 
is shown in Fig. [3] We find the RMS error has the least value 0.12 at p s = 0.3 and the greatest 
value 0.58 at p s — 1. In addition, the RMS errors between p and p man from lfl9ll are nearly the 
same regardless of p s value. 

Table 1 The mean values and standard deviations (<x) of local density and individual velocity in 



measured section, which is gained from the video extraction by the manual procedure IU9I1 



N 


p[l/m] 


p" M "[l/m] 


v m ""[m/s] 


1 






1.24(0.15) 


15 


0.87 


0.77(0.12) 


0.90(0.05) 


20 


1.16 


1.07(0.11) 


0.56(0.05) 


25 


1.45 


1.39(0.12) 


0.34(0.04) 


30 


1.73 


1.71(0.17) 


0.23(0.03) 


34 


1.97 


1.76(0.24) 


0.17(0.03) 



We compare our simulation results with the empirical data (see Table 1 from fl9ll ) and make 
a summary of two groups of results with parameters p s = 1 and p s = 0.3 in the Table 2 and Table 
3 respectively. For Table 2, the residual error Av is the greatest at N = 25 and the smallest at 
N - 15 among five scenes. For Table 3, the residual error Av is the smallest at N = 25 and the 
greatest at N = 15. At low density (N = 15), the difference between the mean velocity of SR 
model and standard LG model is not significant. But the mean gap (d = 1/p) is a little longer 
(Ice// < d < Icell) in SR model than in LG model. At high density (N =30 and 34), the block 
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Table 2 The mean values and standard deviations (cr) of local density and individual velocity in 
measured section by simulation when p s = 1 (standard LG model). The unit of simulation have 
been transformed to physical unit. 



N 


p[l/m] 


p[l/m] 


v[m/s] 


Av[m/s] 


15 


0.87 


1.09(0.00) 


1.24(0.00) 


0.34 


20 


1.16 


1.16(0.00) 


1.24(0.00) 


0.68 


25 


1.45 


1.47(0.05) 


1.04(0.05) 


0.70 


30 


1.73 


1.76(0.06) 


0.87(0.05) 


0.64 


34 


1.97 


1.99(0.02) 


0.67(0.03) 


0.50 



Table 3 The mean values and standard deviations (cr) of local density and individual velocity in 
measured section by simulation when p s = 0.3 (SR model). The unit of simulation have been 
transformed to physical unit. 



N 


p[l/m] 


p[l/m] 


v[m/s] 


Av[m/s] 


15 


0.87 


0.87(0.02) 


1.15(0.03) 


0.25 


20 


1.16 


1.19(0.05) 


0.61(0.03) 


0.05 


25 


1.45 


1.44(0.02) 


0.36(0.01) 


0.02 


30 


1.73 


1.74(0.03) 


0.20(0.01) 


0.03 


34 


1.97 


1.98(0.04) 


0.12(0.01) 


0.05 



states dominate the pedestrian flow for the two models. At medium density (N =20 and 25), the 
moving states of pedestrians in SR and LG model have important difference. For example, all 
pedestrians is in free states in LG model at N = 20. However, a part of pedestrians is in slow 
reaction states or block states in SR model. As an extreme case for p s = 0, no pedestrian can 
move forward in SR model if N > 22. Only leaving two cells in the front can a pedestrian move 
in free state in SR model. 

To validate the distribution of individual velocity and density at fundamental diagram, we 
draw velocity-density diagrams with values of p s =1, 0.3 and 0.1. The local density p, of pedes- 
trian i is the mean value of the momentary density p(f) during the time-slice [?'", t°'" ]. The (v, ,p,) 
distribution of standard LG model in Fig. |4]a) is very different from Fig. IHb), Fig. HJc) and 
Fig. |4|d). The LG model with p s = 1 become a deterministic discrete model and everyone's 
individual velocity at steady state is fixed without fluctuation. For each density, the individual 
velocities distribute along a line. All of the individual velocities of N =15 and 20 are same, which 
reach the free speed (1.24m/s). This is because the critical density p c = 1.27[l/m] (N = 22) in 
our simulation system, after which the stop-and-go waves occur and phase separation is observ- 
able. The discussion about critical density and stop-and-go waves are set up later. In 11911 . the 
density-velocity points for the cycles with N = 15, 20 and 25 distribute in clearly separated areas 
and these areas blend for the cycles with N = 30 and 34 (see Fig. Etd)). However, the values 
of individual velocity distribute broader in our simulation than in 11911 . There are no obvious 
boundaries between the density-velocity points of two different N in Fig.Ub) and Fig.|4|c). A 
possible explanation is that the slow reaction gap with only one cell between the front and back 
pedestrian is not long enough. One observes that the values of individual velocity at low den- 
sity N - \5 and 20 in simulation are higher than that of empirical data. The phenomenon is 
due to the original discreteness of LG model in terms of velocity and position. In our model, a 
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Fig. 4. Dependency between the individual velocity and density of the 100th cycle with N = 15, 
20, 25, 30 and 34. (a) p s = 1 (standard LG model); (b)p s = 0.3 (SR model); (c) p s = 0.1 (SR 
model); (d) Seyfried's experimental data in [ 19]. 



pedestrian either moves one cell with v,(f) = lcell (1.24m/s in physical unit) or stays still with 
v,-(f) = Ocell. There is on mechanism to generate an intermediate value 0cell< v,(f) <lcell about 
the momentary individual velocity. 

In Fig.[2a), the individual velocities do not distribute widely in di - v, diagram but arrange 
in several lines. In Fig. |5jb), a linear relationship with d = 0.46 + 0.63v gives the best fit 
to our simulation data with coefficient of determination R 2 = 0.96. Although the linear trend 
agrees with lfl9ll . the slope 0.63 of fitting line is a little lower than the slope 1.06 in Seyfried's 
experiment. One possible explanation is that the size of cell is square with 0.4m*0.4m. In real 
life, the shape of pedestrian is an approximate ellipse with a typical value of 0.58m*0.33m 131 1. 
Another reason is that due to spatial discretization of LG model, pedestrians cannot be distributed 
uniformly in the grid. 

The stop-and-go waves and phase separation in pedestrian dynamic of standard LG model is 
shown in Fig. [6] With the compression effect of stop-and-go wave, the pedestrian flow is divided 
into two phase: the free phase and jam phase. The common border of the two phases goes 
upstream circularly with propagation velocity v p = -lcell/step (i.e. -1.24m/s in physical unitand 
minus denotes the opposite direction from pedestrian flow). The stop-and-go waves in LG model 
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Fig. 5. Dependency between required length and velocity of the 100th cycle with N = 15, 20, 
25, 30 and 34. (a) p s = 1; (b) p s = 0.3. 



is similar to the stop-and-go waves of highway traffic, which often happens at the upper stream of 
highway bottleneck l32ll . Although the stop-and-go wave and phase separation in the single-file 
pedestrian flow was not reported in lll9ll . the same phenomena were observed and discussed in 
the field experiments of |29|]. A series of larger scale field experiments were arranged by the 
same research team in the wardroom of Bergische Kaserne Diisseldorf at 2006 1E9I1 . The group 
of test persons was composed of soldiers instead of students. The experimental set-up is similar 
to that of lfl9ll except the size of passageway and the number of testers. In addition, similar result 
was reported in simulation by a social force model in POi . We observed in simulation that the 
critical density p c = 1.27[l/m] (corresponding N — 22 in L — 1.73m), after which the stop-and- 
go waves are observable. This coincides with the simulation result of p c = 1.21[l/m] in |[30ll . 
However, one can not observe in standard LG model that some pedestrians leave much larger 
than average gaps in front, which was reported in 1I30I1 . 

The stop-and-go waves and phase separation for SR model is shown in Fig. [7] The global 
stop-and-go wave observed in standard LG model does not exist in SR model. Instead, several 
small density waves go upstream with propagation velocities \v p \ < 1 . For N > 14, the stop- 
and-go waves occur and phase separation is observable. The pedestrian moving in free state 
must leave two empty cells at least in front of itself. Therefore, the critical density p c =43/3 = 
14.3cell in theory. One can observe in SR model that some pedestrians leave much larger than 
average gaps in front, which coincides with [|30h . The pedestrian flow is jostled by the stop-to-go 
waves and forms several jam clusters separated by one or two empty cells. The fluctuation of 
local density and individual velocity indicates the inner randomness of pedestrian movement and 
is in accordance with the observations in IU9I1 and [291. 



4. Conclusions 



In this work, we propose a model called the slow reaction model (SR model) to describe 
the pedestrian's delayed reaction in single-file movement. We simulate and reproduce Seyfried's 
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Fig. 6. (Colored) The propagation (rightward) of stop-and-go waves in the pedestrian flow from 
time step 5,000~5,100 when p s = 1 (standard LG model). The series of red dots denotes the 
trajectory of a pedestrian during this time, (a) iV = 21; (b) N = 22; (c) N = 28; (d) N = 34. 



field experiment at the research centre Jiilich in 2005 111911 and use its empirical data to validate 
our SR model. We compare the simulation result of SR model with that of the standard LG 
model. Our simulation work includes: (a) the fundamental diagram about the relation between 
one-dimensional density and velocity; (b) dependency between the individual required length 
and velocity; (c) stop-and-go waves and phase separation. We test different probability of slow 
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Fig. 7. (Colored) The propagation (rightward) of stop-and-go waves in the pedestrian flow from 
time step 5,000~5, 100 when p s = 0.3 (SR model). The series of red dots denotes the trajectory 
of a pedestrian during this time, (a) N = 15; (b) N = 22; (c) N = 28; (d) N = 34. 



reaction p s in SR model and found the simulation data of p s = 0.3 fit the empirical data IU9I1 best. 
The RMS error of mean velocity between simulation and empirical data is the smallest (s = 0.12) 
at p s = 0.3. The RMS error of SR model p s = 0.3 is smaller than that of standard LG model 
p s = 1 (s = 0.58). In the range of p s = 0.1 ~ 0.3, our fundamental diagram between velocity 
and density by simulation coincides with field experiments in 01911 . The distribution of individual 



10 



velocity in fundamental diagram in SR model agrees with the empirical data Ill9ll better than that 
of standard LG model. In addition, we observe the stop-and-go waves and phase separation 
in pedestrian flow by simulation, which were observed in [30] and 029I1 . We reproduced the 
phenomena of uneven distribution of interspaces by SR model while the standard LG model did 
not implement. The uneven distribution of interspaces was also reproduced by social force model 
in l30ll . The fluctuation of local density and individual velocity indicates the inner randomness of 
pedestrian movement. The SR model can reproduce the evolution of spatio-temporal structures 
of pedestrian flow with higher fidelity to Seyfried's experiments (19] than the standard LG model. 

There are some imperfections for current SR model. The spatial partition is not fine enough 
and the expression of velocity lacks the middle levels. It is hard to describe the individual mo- 
mentary velocity between Om/s to 1.24m/s. This mechanism will make the mean individual ve- 
locity larger at low density and smaller at high density than that of the empirical data. Although 
we introduce the slow reaction mechanism and calculate the mean velocity of several cycles to 
alleviate this problem, this imperfection has not been solved fully until now. Introducing a finer 
discretization of the space for formulating pedestrians' behavior in SR model is is worth further 
research. 
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